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ABSTRACT 

Currently, two competing models are invoked in order to explain the observable prop- 
erties of Anomalous X-ray Pulsars (AXPs). One model assumes that AXP emission is 
powered by a strongly magnetized neutron star - i.e., a magnetar. Other groups have 
postulated that the unusually long spin periods associated with AXPs could, instead, 
be due to accretion. As there are severe observational constraints on any binary accre- 
tion model, fossil disk models have been suggested as a plausible alternative. Here we 
analyze fossil disk models of AXPs in some detail, and point out some of their inherent 
inconsistencies. For example, we find that, unless it has an exceptionally high magnetic 
field strength, a neutron star in a fossil disk cannot be observed as an AXP if the disk 
opacity is dominated by Kramers' law. However, standard alpha-disk models show that 
a Kramers opacity must dominate for the case logB > 12, making it unlikely that a 
fossil disk scenario can successfully produce AXPs. Additionally, we find that in order 
to sufficiently spin down a neutron star in a fossil disk, an unusually efficient propeller 
torque must be used. Such torques are inconsistent with observations of other accreting 
sytems - particularly High Mass X-ray Binaries. Thus, our analysis lends credence to 
the magnetar model of AXPs. 

Subject headings: accretion, accretion disks - pulsars: general - stars: neutron - X-rays: 
stars 



1. INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) are a recently discovered subclass of X-ray pulsators sharing 
distinct properties that are markedly different from that of their high-mass X-ray binary (HMXB) 
and low-mass X-ray binary (LMXB) cousins (Mereghetti &: Stella 1995; Van Paradijs, Taam, &: 
Van den Heuvel 1995). In particular, the six known AXPs have an extremely narrow range of 
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pulse periods, ~ 6 — 12 s. This is remarkable when compared to HMXBs, for example, with spin 
periods in the range 0.069 - ~ 10^ s. Additionally, AXPs have relatively low X-ray luminosities, 
(Lx ~ 10^^—10"^^ erg s~^), no known optical companions and are undergoing stable spindown. They 
also have relatively soft spectra, fit by a combination of power-law and blackbody models. Finally, 
the neutron stars at the heart of AXPs are thought to be young, as evidenced by measured spindown 
ages in the range 10^ — 10^ yr. This assumption is enhanced by the observation of several clear 
supernova remnant associations. A related class of objects (not addressed in this paper) are the 
so-called Soft Gamma Repeaters (SGRs) which share, in their quiescent state, similar spin periods, 
spindown ages, and luminosities. For a recent observational review of AXPs, see Mereghetti et al. 
(2002). 

The major theoretical hurdle in understanding AXPs lies in discovering the source of their 
X-ray emission. Rotational power is clearly not responsible for their power as £^rot = ^ -Z>x 
for typical values. There have been many attempts to model AXPs in recent years but all models 
can, more or less, be grouped into one of two categories. The first class of models assume that the 
source of AXPs are isolated neutron stars with magnetic field strengths in the range 
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Gauss - i.e. "magnctars" (Duncan & Thompson 1992). It was discovered that magnetars could 
account for the timing properties of both AXPs and SGRs by simply invoking standard magneto- 
dipole braking (Thompson & Duncan 1995, 1996). Furthermore, in this model, the source of the 
X-ray luminosity may be explained by rapid magnetic field decay (Thompson &i Duncan 1996; 
Colpi, Geppert, &; Page 2000). The second class of models used to explain the observable features 
of AXPs invoke sonic type of accretion. Such models do not require neutron stars with unusually 
strong magnetic fields; they are assumed to have properties consistent with ordinary radio pulsars, 
B < 10^3 G. In this case, the spindown torque is external and, supposedly, a natural consequence 
of accretion braking. It was quickly shown, however, that traditional accretion models lead to a 
number of problems. Standard Bondi accretion models (Bondi 1952) show that unless AXPs all 
lie near particularly dense regions of the Interstellar Medium (ISM), for example, near molecular 
clouds [e.g. Israel, Mereghetti, & Stella (1994)], it is unlikely that they are directly accreting from 
the ISM. Mereghetti & Stella (1995) suggested that AXPs may be members of very low-mass X-ray 
binaries, (VLMXBs) - a subclass of LMXBs with slightly higher magnetic field strengths, ~ 10^^ 
G. One possible example of a VLMXB is the X-ray pulsar, 4U1626 — 67. Verbunt, Wijers, & 
Burm (1990) had estimated that a low-mass donor star of about O.O2M0 accreting onto an old 
neutron star may be responsible for its X-ray emission and, in fact, an optical identification has 
been made (Chakrabarty 1998) for this system. Although this scenario can not be ruled out for 
AXPs, there are severe observational restrictions that necessarily limit the possibility of detecting 
any potential optical companion they might have (Mereghetti, Israel &: Stella 1998). For example, 
there have been no detected Doppler modulations in the X-ray frequency for any of the AXPs. 
Such modulations are expected, however, even in the case of accretion from a low-mass Helium 
star or a very low-mass white dwarf. Additionally, there seem to be clear observational differences 
between the AXPs and 4U 1626 — 67. For example, the LMXB has a much harder spectrum than 
any of the anomalous X-ray pulsars. Another discrepancy is that while 4U 1626 — 67 has undergone 
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clear intervals of both spinup and spindown (to be expected for an accreting source), this is not 
the case for the AXPs which show persistent spindown only. The first non-binary accretion model 
was suggested by Van Paradijs, Taam, & Van den Heuvel (1995). In this scenario, AXPs are 
descendants of close HMXBs and are observed during the post Thornc-Zytkow/spiral-in phase. In 
this case the source of the X-ray luminosity is a neutron star accreting from a residual disk formed 
from the disrupted companion star. The apparent youth of AXPs is consistent with this model 
as is their association with supernova remnants. Ghosh, Angclini, & White (1997) suggested that 
this model can naturally explain the multi-phase spectra of AXPs. For example, disk accretion can 
account for the power-law phase of the spectra whereas a spherically symmetric flow, emanating 
from the companion's envelope would account for the blackbody component. Of course, this model 
naturally assumes the nascent neutron star can survive the common-envelope/spiral-in phase of its 
evolution. It has been demonstrated, however, that this is not the case (Brown 1995; Chevalier 
1993) as hypercritical accretion forces > IMq onto the Thorne-Zytkow object, crushing it, forming 
a black hole. 

It is unknown whether or not all of the mass ejected from a core-collapse supernova can escape 
the gravitational well of the embedded compact star. With this in mind, it is certainly plausible 
that a significant amount of mass may fall back onto a young neutron star after the supernova (Lin, 
Woosley, & Bodenheimer 1991; Chevalier 1989). Most recently, several groups have independently 
proposed that fossil disks may interact with young neutron stars and that this interaction can result 
in timing signatures similar to what is observed for AXPs (Chatterjee, Hernquist, &; Narayan 2000; 
Alpar 2001; Marsden et al. 2001). In this latest challenge to the magnetar hypothesis it is asserted 
that neutron stars embedded in fossil disks have conventional field strengths, similar to those of 
radio pulsars. Chatterjee, Hernquist, & Narayan (2000) (hereafter CHN) arc the first to present a 
detailed physical model of fossil disk accretion incoporating a time-dependent mass-transfer rate. 
In this model, it is assumed that after a brief transient phase, mass-loss is self-similar, obeying a 
power-law decay. The narrow range of observed spin periods can be explained by assuming the 
inflow quickly becomes an advection-dominated accretion flow (ADAF). Alpar (2001) goes further 
and posits that not only can AXPs (and possibly SGRs) be explained by fossil disk models but 
varying fossil accretion rates can account for all "non-standard" young neutron stars that do not 
appear to be radio pulsars (Kaspi 2000). Marsden et al. (2001) suggest that the environment 
surrounding AXPs (and SGRs) are unusually dense, thereby facilitating the formation of a fossil 
disk. This so-called "pushback" model has been criticized however (Duncan 2002) as the inferred 
densities calculated by Marsden et al. (2001) very strongly depend on the AXP distance, a quantity 
which is not known with accuracy. In any case, the unifying feature of all these fossil disk models 
is that the torque causing the necessary rapid spindown is induced by the propeller-effect (Pringle 
Sz Rees 1972; Illarionov & Sunyaev 1975; Fabian 1975). The propeller mechanism allows a rapidly 
rotating neutron star to prevent the flow of material from accreting to its surface. The resulting 
interaction between the neutron star magnetosphere and the accreting material somehow causes 
the neutron star to lose angular momentum. 
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Here, we show that all fossil disk models carry inherent inconsistencies. We focus our analysis 
on the CHN model, as it is the only current model to provide a detailed physical picture of fossil 
disk accretion but our arguments apply generally. In § 2, we show that the final spin period of 
a neutron star embedded in a fossil disk strongly depends on the disk opacity. We find that, if 
free-free and bound-free transitions dominate over scattering in the disk and B <3.7 x 10^^ Gauss, 
the neutron star can not be observed as an AXP. Employing a simple alpha-disk model (Shakura 
& Sunyaev 1973), we show that, in fact, a Kramers opacity must be assumed if one is to believe 
the fossil disk mass-transport rate of CHN. 

In § 3, we discuss the consequences of the propeller torque model employed by CHN. In 
order to spin the neutron star down to the observed AXP range in a time approximating its 
spindown age, CHN assumed an "efficient" propeller torque, allowing the ejected plasma to be spun 
up to corotation with the pulsar. Although the exact mechanism by which a neutron star loses 
angular momentum by the propeller effect is not well understood, we show that other "standard" 
torque mechanisms, in fact, yield the opposite result. As pointed out by Thompson (2002), if 
one, instead, assumes the plasma is ejected with the specific angular momentum of a particle in 
a Kcplerian orbit at the magnetospheric boundary, one must then assume a dipole field several 
orders of magnitude greater than what is inferred for radio pulsars in order to yield similar results. 
In fact, we will show that such propeller mechanisms fail to produce AXPs. Of course, the reality 
of the magnetospheric interaction is a very complex magnetohydro dynamical problem so it is hard 
to say, a priori, which propeller model is inherently more accurate therefore it is necessary to 
treat all models as competitive and compare the results. In § 4, wc examine the effect of invoking 
various propeller torques to the (hopefully) well-understood HMXB, SMC X-1. We find the CHN 
"efficient" propeller suggests a spin period for SMC X-1 that deviates orders of magnitude from its 
observed period. In § 5, we summarize our findings and discuss the possible implications. 



2. FOSSIL DISK MODELS 

It has long been considered that immediately following the formation of a neutron star in a core- 
collapse supernova, some fraction of the ejected plasma may be unable to escape the gravitational 
well of the star. Following Lin, Woosley, &; Bodenheimer (1991) and Chevalier (1989), CHN suggest 
a fallback mass < O.IMq is not unreasonable. At that point, it is unclear what fraction of this 
fallback mass will eventually form the fossil disk but it is suggested that a large amount can either 
be immediately ejected from or accreted onto the neutron star. Comparing this problem with 
that of the tidal disruption of stars by massive black holes, discussed earlier by Cannizzo, Lee, & 
Goodman (1990), CHN find that fallback matter can circularize into an accretion disk in (roughly) 
a local dynamical time. For their analytical model, CHN assume the dynamical timescale is ~ 1 ms 
but claim that the final outcome is insensitive to the exact numerical value. Furthermore, an initial 
disk mass, Mq = 0.006Mq is assumed. Later wc will discuss the effects of varying the disk mass 
(as well as the initial spin period of the neutron star) but, for now, we note that as the observed 
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spin periods of AXPs are in a very narrow range, it is not likely possible to vary the initial disk 

mass much and find results consistent with the predictions of CHN. The newly born neutron star is 
assumed to have "canonical" properties with a magnetic field strength in the range 10^^ < B < 10^'^ 
Gauss and an initial spin period, Pq ~ 15 ms. Such properties agree with what is observed for the 
Crab pulsar (PSR 0531 — 21) although recent analyses have shown that this may not be the norm 
(Kaspi 2000; Alpar 2001). 

Unlike the case of binary accretion, the fossil disk can not be replenished so accretion is 
necessarily a time-dependent phenomenon. Following Cannizzo, Lee, Sz Goodman (1990), CHN 
suggest that after the brief transient phase, the accretion rate declines self-similar ly. In our analysis, 
we have devised a computer model that can exactly reproduce the analytical model of CHN with 
the added benefit that it is simple to vary any parameter and determine its influence on the final 
outcome. We discuss the results of the model in §2.1 and explore the consequences of varying one 
particular parameter (the opacity) in §2.2. 



2.1. Self-Similar Fossil Disks 

Following the arguments of Cannizzo, Lee, Sz Goodman (1990), CHN suggest that after a 
dynamical time the fossil disk loses mass self-similarly: 




Here, to ~ 10~^ s, is the local dynamical time and F > 1 is a constant that directly depends on 
the disk opacity, k (see below). Not all of the mass lost from the fossil disk will be accreted onto 
the neutron star surface. Thus, in general, the neutron star accretion rate, < M and Mx = 
during the propeller phase. Of course, only the surface accretion rate gives rise to the observed X- 
ray luminosity, = GM^M^/Rx- Normalizing to the total initial disk mass, Mq, it is determined 
from equation (1): 

T- 1\ /Mo'^ 



From Cannizzo, Lee, Sz Goodman (1990), we see the parameter T{k) directly depends on the disk 
opacity. For a standard parameterization, k{p,T) = kqp^T", where p and T are the local density 
and temperature respectively, one finds: 

^ 38 + 18A-4i/ 

32 + 17A-2Z/ ^ ' 

Thus, a standard Kramers opacity {X = l,u = —7/2) yields F = 1.25 whereas if the opacity is 
predominantly due to Thomson scattering (A = = 0), then F = 19/16 = 1.188. In their analytic 
model, CHN suggest F = 7/6 ^ 1.167. 
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The subsequent evolution of the system can be divided into four phases depending on the 
relative strengths of several parameters. The magnetospheric boundary of a neutron star is often 
defined as the point where ram pressure of infalling matter is balanced by the pulsar's magnetic 
dipole pressure [see Francischelli, Wijers, & Brown (2002) and references therein]. Hence, the flow 
of accreted matter is governed by the dipole field for any r < i^m- Here i^m is the magnetospheric 
radius which, for a neutron star with = 10 km and = IAMq, is determined to be: 



In equation (4), MEdd ~ 9.46 x 10^^ g is the (hydrogen) Eddington accretion rate, Bu = 
Gauss, and the factor of 1/2 comes from the assumption of a disk geometry. Regardless of whether 
the flow starts out in a disk or has a spherical geometry, gravitationally captured plasma can only 
reach the neutron star's magnetosphere if it can fall along a closed field line. This is only permitted 
within the neutron star's light cylinder, defined such that Ric = c/il = cP jlix. Thus, plasma can 
interact with the neutron star's magnetosphere if and only if i?m < i?ic. If this is not the case, the 
neutron star evolves independently of the fossil disk and spins down by magnetodipole braking (i.e. 
emitter phase). Once the neutron star spins down sufficiently such that matter can couple to its 
magnetic field, the interaction's effect on the overall spin evolution depends on the balance between 
centrifugal and gravitational accelerations. If the pulsar is initially rotating too fast, the neutron 
star will eject infalling plasma, propelling it away tangcntially, while simultaneously losing angular 
momentum in the process. This "propeller" effect (Pringle & Rees 1972; Illarionov &: Sunyaev 1975) 
can be parameterized by a fastness parameter Ug = r2/r2K(-Rm)) where ClK{Rm) = {GM^/R'^)^/^ is 
the Keplerian angular velocity at the magnetospheric boundary. If the fastness parameter greatly 
exceeds unity, the propeller mechanism is initiated. For both the propeller and emitter phases, 
no material can accrete to the neutron star surface and oc ~ 0. As discussed in the 
introduction, the way in which a propelling neutron star loses angular momentum is not well 
understood. However, it is safe to assume that some kind of propeller torque (see equation [6], §3) 
acts on the neutron star, spinning it down with time, and reducing the fastness parameter. When 
Us < 1, surface accretion can take place and, provided the transfer rate is sub-Eddington, M = M^. 
Contrary to what is expected in an X-ray binary, however, M steadily decreases in a fossil disk 
and an extended accretion phase never really occurs. Instead, a quasi-equilibrium period, which 
CHN dub a "tracking" phase, occurs, and it is during this period that the neutron star can be 
observed as an X-ray pulsar. Furthermore, they suggest that the tracking phase must be short 
lived as the accretion flow quickly becomes advection-dominated (ADAF). It is proposed that an 
ADAF transition occurs when the accretion luminosity falls to lO^^Lsdd ~ 1-8 x 10^^ erg s~'^ 
(Chatterjee, Hernquist, & Narayan 2000; Narayan & Yi 1995). At this point, it is thought that 
most of the captured matter is ejected prior to reaching the neutron star's surface and, once again, 
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Mx ~ 0. Prom equation (1), the ADAF transition time is given by 



( 



lOOMo 

AfEdd 



) 



1/r 



iADAF — td 



(5) 



Typically, tADAF ~ 20, 000 — 40, 000 yr. In sum, the neutron star can only be observed as a bright 
X-ray source during the tracking phase of its evolution, ttrack 2tADAF- 

To determine the spin evolution of the neutron star + fossil disk system, a spinup/spindown 
torque is needed. For < Rm < Ric, CHN propose the following (Menou et al. 1999): 



Assuming an arbitrary value for F, equations (1) and (6) can be combined to yield an analytic 
formula for 0,(t), in terms of incomplete gamma functions. For the particular case F = 7/6, CHN 
find a solution using exponential integrals (see their eq. 4). Instead of using their analytic method, 
however, we solved the above equations numerically, using an Eulerian integration scheme. The 
results of our analysis (with F = 7/6) is shown in Figure 1. Here we plot the spin evolution of 
the neutron star as a function of magnetic field strength until the approximate onset of the ADAF 
phase, 2tADAF ~ 38, 000 yr. Initial conditions Pq = 15 ms and Mq = 0.006Mq are assumed. We 
see that a high magnetic field strength is needed in order for the neutron star to be observed as an 
AXP. For the above initial conditions, neutron stars with B12 < 3.9 do not spin down efficiently 
and, consequently, are observed as radio pulsars. On the other hand, for B12 > 3.9, the neutron 
star, after a brief emitter phase, undergoes a rapid propeller spindown and enters the tracking 
phase of its evolution at ttrack ~ 10^ yr- Our numerical analysis agrees with the analytic model of 
CHN (see their figure 1). 

The initial spin period of neutron stars is largely unknown and the long-accepted paradigm 

that assumed the Crab pulsar is the prototypical young neutron star is now being challenged (Kaspi 
2000). Thus, it is important to investigate how the previous results depend on Pq. In fact, we find 
that the duration of the AXP phase does not strongly depend on the initial period. For example, 
assuming Pq = 150 ms had little effect on the final evolutionary state of the system (pre- ADAF). 
For a particular case, consider a neutron star with B12 = 7.5 and Pq = 150 ms. Our analysis has 
shown that it is able to spin down via the propeller effect to the AXP range (10.0 s) in a time 
^ADAF = 4 X 10^ yr. The only deviation from what is seen in figure 1 is in the early evolution. In this 
case, we found a short-lived accretion phase can occur ~ 100 yr after birth. This is compensated 
by a somewhat steeper propeller cycle and a slightly shorter tracking phase. The overall evolution, 
however, is qualitatively the same as the 15 ms case and, as this case is representative of the general 
trend, we therefore conclude that the CHN model can accomodate a wide range of initial periods. 

Another free parameter to be considered is the initial disk mass. Even assuming a fallback mass 
< O.IM0, it is unclear how much matter will circularize into an accretion disk. In the initial period 
following the supernova, much of this material may be accreted onto the neutron star, possibly 



j = in = 2MR'^nKiRm) 1 




(6) 
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hypercritically. Additionally, as detailed modelling probably depends on the progenitor history, we 
make only a few qualitative remarks about assigning possible values to the disk mass. As CHN 
have shown, an AXP can be formed even if Mq <^ O.IMq. However, there still must be a sufficient 
amount of disk matter in order to penetrate the neutron star's light cylinder. From equation (4), we 
see the magnetospheric radius increases as the accretion rate falls and if M falls below some critical 
value, matter will be unable to couple to the closed field lines and no (extended) propeller phase will 
take place. Our numerical simulations indicate that Mq > 7.5 x lO^^M© is needed for the neutron 
star to be observed as an AXP. A neutron star embedded in a fossil disk with significantly less mass 
than this, rarely leaves the emitter phase of its evolution and, hence, would be observed as a radio 
pulsar. The exact critical value strongly depends on the neutron star's magnetic field strength; 
for Bi2 > 7.5, even more disk matter is needed. The most important observational constraint in 
determining the initial disk mass is the narrow distribution of AXP spin periods. Deviations far 
from the CHN sample value of O.OOGMq, although often producing neutron stars with a tracking 
period, often yield a wider range of spin periods than what is currently observed. Thus, for our 
model, Mq = O.OO6M0 was consistently used. 



2.2. The Effect of Varying Disk Opacity 

One quantity that deserves special attention is the disk opacity, parameterized by T(k) in 
equation (3). We now show that if free- free and bound- free transitions dominate over electron 
scattering within the disk (i.e., if a Kramers opacity is assumed), then a canonical neutron star will 
never reach a tracking phase and, consequently, will not be observed as an AXP. We recall from 
equation (3) that Fes = 19/16 = 1.1875 and Fxr = 1-25 so it is not immediately obvious that such 
numerically similar parameters could lead to completely different evolutionary scenarios. Figure 2, 
however, illustrates a strong correlation between F and the overall spin evolution of the neutron 
star. Here we plot the final spin period of a neutron star embedded in a fossil disk as a function 
of F(k) for various magnetic field strengths. An efficient propeller torque given by equation (6) is 
assumed. As we have seen in section § 2.1, the system's evolution also depends on the neutron star's 
magnetic field. In particular, for a given disk opacity, the field strength determines the extent of the 
propeller phase and, ultimately, whether or not a tracking phase can occur at all. This dependence 
is shown in figure 2 by comparing opacity-period evolution curves for B12 = 2.5, 5.0, 7.5, 10.0 and 
40.0. Although the shape and qualitative behavior are the same for the various field strengths, the 
overall range of final periods differs greatly. Generally, each curve shows essentially no dependence 
on opacity below some cutoff value. As F increases past a value, Fc, however, there is a sharp 
drop in final spin period followed by a more gradual decline. These curves illustrate the fact that, 
depending on the strength of the field, a neutron star embedded in a disk with opacity parameter 
greater than some Fc, will not be able to spin down to the AXP range and, consequently, it will 
end its observable life as a radio pulsar. The decline in final spin period with increasing F, can 
be interpreted physically by noting from equation (1) that the overall mass-transport rate steeply 
decreases with time by a factor oc and for F > Fc, disk matter is depleted too quickly to 
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sustain an extended propeller phase. A secondary cause for the decline is the ADAF transition, 
and subsequent X-ray shut-off, must necessarily occur at earlier times for increased T. For T = 7/6, 
our results agree with the conclusions of CHN. As mentioned previously, in this case, the minimum 
field strength needed to support a tracking phase ~ 3.9 x 10^^ G. Thus, in figure 2, Fchn < Tc for 
all fields except B12 = 2.5. If the fossil disk were to be dominated by electron scattering, we find 
that a stronger field is needed in order for the neutron star to enter a tracking period before the 
ADAF phase of its cvohition. In this case, the minimum field required is ~ 6.6 X 10^2 Gauss. For 
B12 = 5.0, Fes > Fc and for B12 = 7.5, Fes Fc. If the disk is dominated by a Kramers opacity, 
however, unusually strong magnetic fields are necessary in order for the neutron star to be observed 
as an AXP. In fact, Fxr 3> Fc for all field strengths below ~ 3 x 10^^ Gauss. All neutron stars 
with magnetic fields below Bxr ~ 3.7 x 10^^ Gauss and otherwise similar initial conditions will 
remain a radio pulsar throughout its observable lifetime. Of course, the strength of this conclusion 
directly depends on the exact point at which the ADAF transition commences. As this quantity is 
not known with accuracy, the true value of Bj^r may be either higher or lower. 

We have seen that the fossil disk opacity is a key parameter in determining the ultimate 
evolution of the neutron star. In particular, we have shown that if a Kramers opacity dominates, 
it is unlikely (though not impossible) for an AXP to be observed whereas if electron scattering is 
more important, it is more likely that a "canonical" neutron star can become an AXP. In fact, 
we now show that, for B12 > 1, a Kramers opacity must be assumed. In order to calculate 
the opacity we adopt the so-called alpha-prescription of Shakura & Sunyaev (1973), implicitly 
assuming that a thin-disk structure is valid for fossil disks. Such an assumption seems valid as it is 
probably reasonable to expect the fossil disk is approximately thin in the regime, t^yn t < tADAF- 
Furthermore, we assume the (kinematic) turbulent viscosity, u, can be parameterized according to 
the standard relation, v = acgh. Here, Cg is the local (isothermal) sound speed and h represents the 
disk's vertical scale-height. It has been argued on physical grounds that the viscosity parameter, 
a < 1 (Shakura &; Sunyaev 1973). The thin-disk assumption is parameterized by the condition 
h/r ~ Cs/v^ <C 1, where r is the radial measure of the disk and v^{r) = r^Kir) = \l GM^/ r 
represents the azimuthal flow within the disk. The general procedure for determining the time- 
evolution of a viscous disk undergoing mass-transfer variations has been extensively studied. Bath 
k Pringle (1981), for example, have shown that the time-dependent conservation equations lead 
to a non-linear diffusion equation for the surface density, S(r, t). It has been shown, however, 
(Prank, King &; Raine 1992), that if at each point in time, and for a given M{t), the timescale 
in which matter diffuses through the disk greatly exceeds the dynamical timescale, ~ r/v^{r), the 
disk may be considered to be approximately steady-state. Except at early times, this condition is 
generally valid for the fossil disk we are considering and, thus, in our analysis, we assume steady- 
state solutions may be used for each value of t. This approach is roughly equivalent to numerically 
integrating the diffusion equation of Bath &; Pringle (1981). 

We first assume an opacity dominated by bound-free transitions, Kxr — 4x 10^^ Z{l+X)pT~'^/'^ 
cm^ g~^, (Schwarzchild, M. 1958) where X and Z < 1 are the hydrogen and heavy element 
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mass fractions, respectively. Following the analysis of Prank, King & Raine (1992), we can then 
estimate the region within the disk where a Kramers-type opacity may be dominant. 



«Kr(r,t) = nZ{l + X)]y^M-'/\t)r%'f-\r) cm^g" 



(7) 



Here M17 = M/IO"*^^ g s^^, and rio = cm as usual and /(r) = [1 — (ri/r)^/^]-*^/^. Fortunately, 

the result is independent of a and, thus, does not depend on the poorly understood viscosity 
mechanism within the disk. The quantity, represents the inner edge of the disk and for neutron 
stars, n = Rui, given by equation (4). In our analysis, we ignore the perturbations brought about 
by edge effects and concentrate on the region, r S> n, /(r) « 1. At temperatures > 10^ K, the other 
dominant source of opacity is electron scattering, k^s = 0.2(1 + X) cm^ g~^. Prom equation (7), 
we can estimate the region of the fossil disk where kkt ^ Hes- 



rKr{t) > (3 X 



106) 



1 + 



M^{^{t) cm 



(8) 



A IAMq neutron star is assumed throughout our analysis. Prom equation (4), we estimate the 
magnetospheric boundary of the neutron star, R^, as the point at which the disk truncates. We 
use this result along with equation (8) to determine the conditions necessary such that electron 
scattering may be ignored. I.e., we calculate the necessary conditions such that Rmit) > ?'Kr(0- 



Mnit) < 50 



B 



3/5 
12 



(9) 



Pinally, combining the above relation with the mass-transfer equation, (1), we can estimate the 
time after which electron scattering may be completely ignored within the disk. Prom (1), (2), 
and (9), we find 



tKr > 20 



B 



-12/25 



12 



(10) 



Por the CHN AXP-cutoff B12 ~ 4, this corresponds to a value of ~ 10 yr. With B12 ~ 37, 
tKr ~ 3.5 yr. In any case, we see that, regardless of the magnetic field strength, the fossil disk 
model of CHN must be Kramers-dominated throughout its evolution and P = 1.25 is the proper 
choice of exponents in equation (1). Coupling this with the preceding analysis and figure 2, we 
have now shown that, unless the neutron star has an extraordinarily high magnetic field, in excess 
of ~ 3.7 X 10^"^ Gauss, the CHN model fails to produce an anomalous X-ray pulsar. 



3. PROPELLER TORQUES 

In its most general form, the exact manner in which angular momentum is transferred between 
an accretion disk and a neutron star is a complex magnetohydrodynamical problem to which no 
simple analytic solution has been determined. More than two decades of numerical analyses have 
resulted in only minor modifications to the original model of Ghosh &; Lamb (1979) [but see Wang 
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(1987)]. It has been determined that when modelhng disk accretion, it is important to include the 
effects of magnetic torques within the disk to the overall spin rate. Ghosh & Lamb were able to 
show that for slow rotators [Q/0,}:q{R„^) = Ug <^ 1], magnetic coupling may enhance spin-up torques 
by as much as 40%. For < 1, the opposite is true, and magnetic effects might actually oppose 
the spin-up. However, as we have seen, due to a time-dependent mass-transfer rate, a neutron 
star in a fossil disk never really has an extensive accretion period nor can it reach true equilibrium. 
Unfortunately, unlike the case for accretion, angular-momentum transfer during the propeller phase 
is not well modelled at all. At present, no ab initio theory exists to compute the torque from an 
accretion disk on a magnetized star. However, over time, several approximate schemes have been 
introduced, generally based on basic physical principles such as conservation laws (Pringle & Rees 
1972; Illarionov & Sunyaev 1975; Fabian 1975). Additionally, Menou et al. (1999) have introduced 
the "efficient" propeller torque, given by equation (6). As several competing models exist, the best 
we can do is compare the models with each other and, most importantly, with observations. With 
this in mind, it is important to emphasize that while we can not really claim one model to be 
intrinsically "better" than another, we can reject a particular model if observations force our hand. 

For a strong propeller, cjg 3> 1, and equation (6) can be written in the form: 

A^CHN = m = -2MRlp. (11) 

Thus, we see the CHN model assumes a rapidly rotating neutron star propels matter from the 
magnetospheric boundary with (twice) the angular velocity of the neutron star, thereby spinning 
it up to corotation. Historically, other propeller models that have been used are not quite so 
efficient. In fact, it has generally been assumed that the propelled matter is, instead, ejected from 
the magnetospheric boundary with the specific angular momentum of a particle in an escaping 
parabolic orbit. 

Nj = m = -MR^v,sc{Rm) = -V2MRlnK{Rm) (12) 

The J-subscript in equation (12) indicates we have used angular momentum methods to estimate 
the propeller torque. An alternative approach is to employ energy conservation principles. Fabian 
(1975) has suggested that, over time, the rotational kinetic energy of the neutron star will be 
transmitted through shocks to the wind plasma falling near the magnetospheric boundary. Conse- 
quently, the infalling gas will heat up and be dispersed when it attains escape velocity. Thus, we 
find E = lOri = — iMt;gg(,(i?in) = —GMM^/R^. The propeller torque is therefore expressed as 

N^ = m = -^i^M£^ (13) 

where uj^ is the fastness parameter. Upon comparison of equations (11) - (13), we find the relation- 
ship between the various proposed propeller torques: 

A^CHN = V2a;sA^j = 2uIn^ (14) 

In the evolution models just discussed here, we found that the fastness parameter varies in the 
range ~ 20 — 75. It is therefore not very surprising that the choice of propeller torques has a great 
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deal of influence on the overall evolution of the neutron star /fossil disk system. This is illustrated in 
figure 3. Here, we plot a period-time relationship for the various propeller models, keeping all other 
initial parameters fixed. We have assumed, as usual, an initial fossil disk mass of Mq = O.OOGMq 
and the initial field strength and spin period are 7.5 x 10^^ Gauss and 15 ms, respectively. In order 
to make our argument independent of the results of the preceding section, we have kept the opacity 
parameter fixed at T{k) = 7/6, the original CHN value. Clearly, employing the efficient propeller 
torque of Chatterjee, Hernquist, & Narayan (2000) makes all the difference in determining the 
overall evolution. When the torque given by equation (11) is used, the neutron star, after a sharp 
propeller cycle reaches a tracking phase at ~ 1.5 x 10^ yr. It is then visible as a bright AXP until 
2iADAF ~ 3.8 X 10^ yr, at which point its final spin period is 13.5 sec. For both the energy and 
angular momentum propellers, the neutron star never reaches equilibrium and despite much higher 
fastness parameters (~ 70), the neutron star's final period never exceeds ~ 0.2 s. Note the early 
evolution for the energy and angular momentum propellers deviate slightly as the energy propller, 
for these initial conditions, lacks an early emitter phase, and begins propelling matter near t ~ 0. 
At later times, the evolution traces the angular momentum propeller case exactly. Of course, if 
given sufficient time, we would see the angular momentum propeller force much higher spin periods 
but only at times t S> 2tADAF- Thus, neither the angular momentum nor the energy propellers 
yield AXPs but, instead, ordinary radio pulsars. 

We recall, from § 2.2, that the choice of r(K) strongly determines the evolution of a neutron 
star in a fossil disk provided the CHN torque is employed. In fact, we found that exceptional 
magnetic fields are required if one were to assume a physically realistic Kramers opacity. The exact 
value of r is much less important to the overall evolution for the other propeller torques, however. 
Assuming a disk dominated by electron scattering, for example, does not change the final spin 
period for either the angular momentum or energy propellers. For a Kramers disk, (F = 1.25), we 
also see no early evolutionary changes in spin period. However, in this case, the onset of ADAF 
flow occurs much earlier (tADAF — 2650 yr by equation [5]). For a Kramers disk, therefore, we find 
the final spin period for all three propeller torques is the same: about 0.13 s. As pointed out by 
Thompson (2002), increasing the magnetic field strength by an order of magnitude does allow the 
less efficient propller torques to spin down the neutron star to longer periods, but we find that, 
regardless of the strength of B, the J-propeller and E-propeller can never produce an anomalous 
X-ray pulsar. For example, in order to spin down a neutron star in a fossil disk by the angular 
momentum propeller torque to a final period of ~ 5 s, it is necessary to increase the initial field 
strength to = 9.8 x 10^^ G. In fact, even at such high fields, the neutron can never reach equilibrium. 
Instead, as ujs ^ 24 (and is slowly decreasing) at t ~ 2tA_BAF, there will be no tracking period at 
all. Borrowing from the terminology of Chatterjee, Hernquist, &: Narayan (2000), we dub such a 
neutron star a "dim propeller." Spun down into the graveyard (Pf ~ 5.3 s), the neutron star is a 
weak emitter in both radio and X-rays. This trend is even more obvious for the case of the energy 
propeller. Here, in order to spin the neutron star down to a final period of 5.9 s, a magnetic field 
> 8.6 X 10^^ G is needed. Even ignoring the fact that such a field strength is already well into 
the so-called "magnetar" range, we find this neutron star, too, is destined to be a dim propeller 
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[u;s(2tADAF) ^ 140]. Thus, we see the results of Chatterjee, Hernquist, & Naxayan (2000), namely 
that AXPs can be produced from ordinary radio pulsars, hinges upon the condition of an efficient 
propeller torque. If one can show that a canonical neutron star can not spin propelled matter at the 
magnetospheric boundary to corotation, then it is a necessary consequence that fossil disk models, 
as outlined by Chatterjee, Hernquist, & Narayan (2000), cannot produce anomalous X-ray pulsars. 

4. A COMPARISON OF VARIOUS PROPELLER TORQUES FOR SMC X-1 

Initially discovered during a rocket flight (Price et al. 1971) and identified by Uhuru as a 
discrete X-ray source in the Small Magellanic Cloud (Leong et al. 1971), SMC X-1 is the most 
luminous and one of the most extensively studied of the X-ray pulsars. Its binary nature, implied 
by periodic X-ray eclipses (Schrcier et al. 1972) was confirmed by the identification of an optical 
companion: BOI supergiant, Sk 160 (Webster et al. 1972; Liller 1973). The 3.89 day orbital period 
is decaying at a rate -Porb/-Porb = (—3.36 ± 0.02) x 10~® yr~^ (Levine et al. 1993), presumably due 
to tidal interactions between the orbital motion and the primary's rotation. The orbital elements 
of the system have been measured with the neutron star's mass estimated at M2 = 0.8 — l.SM© and 
the primary mass and radius given by Mi = 19Mq and -Ri = I8-R0 respective (Primini, Rappaport, 
Sz Joss 1977). For consistency, we will continue to assume M2 = 1.4Mq and R2 = 10® cm. The orbit 
is circular [e < 0.00004, (Levine et al. 1993)] and Kepler's third law suggests an orbital separation 
a ~ 30Rq . With an 0.71 s pulse period (Lucke et al. 1976), SMC X-1 is the only X-ray pulsar 
undergoing stable spinup at a rate P = —1.2 x 10~^^ s s~^ (Kunz et al. 1993; Kahabka & Li 1999). 
Spinup implies the neutron star is accreting and, indeed, there is evidence of an accretion disk 
around SMC X-1 (Van Paradijs 8z Zuiderwijk 1977; Tjemkcs, Zuidcrwijk, &; Van Paradijs 1986). 
Its luminosity, L^, has been measured to vary from ~ 10'^^ erg s^-*^ in the low-intensity state to the 
extremely high value of 5 x 10^^ erg s~^ ~ 5LEdd in the high-intensity state. The presence of an 
accretion disk, excess X-ray luminosity and persistent spinup, together suggest mass-transfer may 
be due to atmospheric Roche-lobe overflow from the massive primary (Savonije 1979). Recently, an 
X-ray burst was discovered from SMC X-1 (Angelini, Stella, &; White 1991), which was suggested 
to be due to instabilities in the accretion flow, i.e., a type II burst (Li & Van den Heuvel 1997). 
If, indeed, SMC X-1 is a member of the class of "bursting pulsars," like the recently discovered 
transient X-ray pulsar, GRO J1744-28, then a diminished pulsar magnetic field might be the cause 
of such bursts. In fact, Li & Van den Heuvel (1997) have estimated a magnetic moment for SMC 
X-1 to be somewhat lower than normal, /j, = BE^ ~ 10^^ G cm^. This is consistent with claims 
that magnetic field decay is connected with binary accretion [see Francischelli, Wijers, &; Brown 
(2002) for a discussion]. 

In order to test the various propeller torques of the preceding section, we apply a modified 
accretion model to the SMC X-l/Sk 160 system. Since the giant phase is only about ~ 10% of the 
main sequence lifetime, it is important to consider the mass transfer mechanism for the progenitor 
system as well. Most likely, the progenitor to the supergiant primary is an 09 V star, with ZAMS 
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mass and radius, ~ 20Mq and ~ IORq. Its main sequence lifetime can be estimated by the mass- 
luminosity relation, L oc M^-^, such that Tms ~ 10^° yr {M/Mq)~'^-^ = 5.6 x 10^ yr. We assume 
the progenitor lies within its Roche lobe throughout the hydgrogen core-burning phase so mass 
transfer is, then, most likely due to a steady stellar wind, of order ~ 10~'^ Mq yr~^ (Kudritzki &: 
Puis 2000). Thus, the overall evolution of the system has two distinct parts: In the first part, an 
09 V star transfers mass, by steady stellar wind, to a young neutron star with initial spin and 
field strength, ~ 15 ms and ~ 10^^ G. During this phase of the evolution, mass transfer is highly 
non-conservative as only a fraction of the mass emitted by the primary can be captured by the 
neutron star, M2 = —fc^i- The non-conservative phase lasts the entire main-sequence lifetime of 
the primary. In the second part of the evolution, a BOI supergiant conservatively transfers mass via 
Roche-lobe overflow (while non-conservative wind-transfer continues) for a time ~ O.Itms ~ 6 x 10^ 
yr. Accretion-induced field decay is assumed for both parts of the evolution. The evolution code 
for such a model, while not much different from the original code for the fossil disk model described 
above, is almost identical to the analysis described in Francischelli, Wijers, & Brown (2002). Thus, 
for details of the input physics assumed, we refer the interested reader there. For the sake of brevity, 
here we simply show the results of our analysis in figure 4. We plot the neutron star spin period 
as a function of time for all three propeller models. Of course, we have no way of estimating how 
long SMC X-1 has been undergoing shell burning so the most we can do is compare the known spin 
period of 0.71 s to the range of predicted values calculated during the supergiant phase. 

If an energy propeller torque is assumed we see that, after a brief emitter phase, the neutron 
star in SMC X-1 undergoes an extended propeller cycle. Weak propcUcr-torquing allows the neutron 
star to gradually spin down until it reaches a maximum period of about 30 s at t ~ 3.1 Myr. At this 
point, accretion spinup commences and continues until a renewed propellcr-spindown at t ~ 4.4 
Myr. This leads to a more gradual accretion cycle and, at the onset of hydrogen-shell burning, at 
5.6 Myr, the neutron star is still in the graveyard at P = 5.3 s. Overall, during the main-sequence 
evolution, the radiation-driven wind can only transfer 3.6 x 10~^Mq to the neutron star surface. 
As soon as the primary overflows its Roche-lobe, however, heavy accretion ensues, recycling the 
neutron star back into observability, spinning it up and bringing down its magnetic field. After a 
total evolution time of 6.2 Myr, Pf = 0.99 s. Employing the accretion-induced field-decay model 
of Francischelli, Wijers, k, Brown (2002), the final magnetic field is calculated to be Pf 4 x 10^° 
G. For the case of an angular momentum (J) propeller, propeller-spindown is much more rapid 
and occurs much earlier in the system's evolution. As a result, the accretion phase also occurs 
much earlier and the pulsar spins back up to ~ 2.2 s at the end of the core-hydrogen burning 
phase. Once again, a rapid accretion cycle occurs during Roche-lobe overflow and the flnal spin 
period is 0.52 s. As the total accretion time is longer, we flnd the final magnetic field is necessarily 
lower for the angular momentum propeller and = 1.3 x 10^° G. Finally, we look at the case of 
the efficient, CHN propeller torque. Here, propeller spindown is even more rapid than before and 
the neutron star reaches P = 31.3 s in only ~ 3.3 x 10^ yr. Unlike what was observed for the 
other two propeller models, here there is no subsequent recycling and the neutron star remains, 
roughly, in spin equilibrium for the duration of the main-sequence evolution. Consequently, no 
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mass is accreted to the neutron star surface during this period. As soon as the primary begins to 
overflow its Roche lobe, however, a new propeller cycle spins the neutron star down further into 
the graveyard, reaching Pmax ~ 100 s. Finally accretion is allowed, but the pulsar can never be 
sufficiently recycled back into observability and Pf = 15.9 s. Altogether, only 9.7 x 10~^Mq is 
accreted to the neutron star surface and, as a consequence of such weak accretion, = 1.2 x 10^^ 
G. 

In conclusion, we see that both the energy and angular momentum propellers, for the given 
initial conditions, yield reasonable spin histories for our model of SMC X-1. The magnetic field 
evolution is not as accurate, however, as the heuristic field-decay model of Prancischelli, Wijers, & 
Brown (2002) , while yielding results consistent with observations of LMXBs and relativistic binary 
pulsars, seems to deviate by a factor or ~ 2.5 — 10 from the magnetic moment calculations of Li & 
Van den Heuvel (1997). Although for the purpose of this paper, wc arc only interested in the spin 
evolution of neutron stars for the various propeller torques, the larger subject of pulsar magnetic 
field decay in general remains an open issue. The most important result of our analysis, is that 
we have definitively shown that the CHN propeller torque model is completely inconsistent with 
observations of SMC X-1. In particular, we have shown that such a propeller torque yields a neutron 
star with a spin period somewhere in the range 15.9 — 100 s, orders of magnitude different from the 
observed period. On this basis, we reject the efficient propeller torque model and determine it is 
extremely unlikely that a neutron star can propel matter to corotation. Since we have also shown 
that less efficient (yet more physically plausible) propeller torques can not spin neutron stars in 
fossil disks down into the AXP range, regardless of their magnetic field strength, we must reject 
fossil disk models of anomalous X-ray pulsars. 



5. DISCUSSION 

Accretion from a fossil disk is the latest scenario that has been suggested as a plausible alterna- 
tive to ultramagnetized neutron stars, or magnetars, as a way of explaining the timing signatures 
of anomalous X-ray pulsars. Although several groups have independently postulated fossil disk 
accretion models, (Chatterjee, Hernquist, Sz Narayan 2000; Alpar 2001; Marsden et al. 2001) only 
Chatterjee, Hernquist, & Narayan (2000), CHN, have provided a detailed, time-dependent mass- 
transfer mechanism and, as a result, our analysis focusses on their model. In fact, we have shown 
that the CHN model contains inherent inconsistencies and, ultimately, when measured against ob- 
servations and realistic phyiscal assumptions, fails. In section 2.2, for example, we have shown that 
a realistic estimate of the disk opacity forces us to either assume unusually strong neutron star 
magnetic fields, or abandon the CHN model altogether. Using the alpha-disk model of Shakura 
Sz Sunyaev (1973), we have shown that a Kramers opacity must be assumed for any fossil disk 
with a mass-transfer rate given by equation (1). This is in contrast with the analysis of CHN, 
where a non-physical opacity parameter, T{k) = 7/6, was assumed. We can interpret this result as 
being due to the fact that the magnetospheric boundary increases with increasing field strength. If 
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the magnetic field, therefore, exceeds some minimum value, the inner edge of the fossil accretion 
disk necessarily truncates far enough away from the neutron star such that temperatures are low 
enough to ignore electron scattering. Furthermore, we have shown that a neutron star in a Kramers 
disk (Fxr = 1-25), can not spin down to the AXP range fast enough for B12 ~ 1 — 30. We have 
estimated the minimum pulsar field needed to give results consistent with CHN to be ~ 3.7 x 10^^ 
G. Although seemingly unusual, there have been recent discoveries of radio pulsars with inferred 
field strengths in this range. Camilo ct al. (2000) have reported the discovery of J1119 — 6127, a 
1670-ycar old pulsar apparently centered in a previously uncatalogued supernova remnant, with an 
inferred field strength B = A.l x 10^^ G. Additionally, the same group has reported the discovery 
of PSR J1814 - 1744 with Tc = 85 kyr and B = 5.5 x 10^^ G. Note that such field strengths are 
quite close to the so-called magnetar critical field. Be = rn^c? /eh = 4.4 x 10^^ G. It has been 
suggested that, dX B > Be, QED processes such as photon splitting may inhibit pair-production 
cascades near the magnetic poles and, neutron stars with field strengths in this range should be 
radio-quiet (Baring & Harding 1998). In fact, both J1119 - 6127 and J1814 - 1744 seem to be 
ordinary radio pulsars and have no discernible X-ray emission. Complicating matters further, we 
note that 1814-1744 has spin properties quite similar to the AXP, IE 2259 -|- 586. Thus, it seems 
that the quantum demarcation line outlined above may be a bit fuzzy and it is not clear whether 
or not a neutron star in a fossil disk with B12 > 3.7 can be observed as a radio pulsar or AXP. 
In any case, the thesis of CHN, namely "ordinary" neutron stars can produce AXPs, seems a bit 
strained at this point. 

Another apparent inconsistency with the CHN model is their assumed propeller torque, given 
by equation (6). We have seen that, for rapid rotators, this model assumes that the neutron star 
forces propelled material to be ejected from the magnetospheric boundary with the same angular 
velocity as the star itself - i.e. it must be spun up to corotation. After noting that such a 
propeller torque is somewhat inconsistent with historical precedent, we realize there is insufficient 
observational data on the propeller effect to immediately rule it out. The exact nature of the torque 
is, however, crucial in determining the overall evolution of neutron stars in fossil disks, as was 
pointed out in section 3. For example, if one assumes material is propelled from the magnetosphcre 
with the angular velocity of a Kepler particle then an AXP will not be observed. Instead, the 
neutron star will either be observed as a radio pulsar, or in the case of an extremely high magnetic 
field {B > 10^^ — 10^^ G), a dim propeller. As a test, we applied the CHN propeller torque 
to the well-studied HMXB, SMC X-1 and compared it to known models of angular momentum 
transfer, namely the energy and angular momentum propeller torques. We foTind that, during 
the supergiant phase, the CHN model predicts a range of spin periods for SMC X-1 ~ 2 orders 
of magnitude different from the known period. It seems that the propeller torque employed by 
Chatterjee, Hernquist, & Narayan (2000) is simply too efficient to be consistent with observations. 
A more robust comparison of propeller torque models would be illuminating. For example, one 
could examine the evolutionary history of all HMXBs with known spin periods. As there is a wide 
range of spin periods observed for these systems (especially the Be HMXBs) , it is unknown whether 
or not one propeller model can account for such variety. In any case, our preliminary work seems 
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to rule out the CHN propeller torque as physically unrealistic and, consequently, strikes a blow 
against fossil disk models in general. 

Ultimately, only observations will definitively determine the true nature of anomalous X-ray 
pulsars. In particular, optical and/or infrared studies of the X-ray sources could help solve the 
mystery of AXPs once and for all. A direct detection of disk emission, for example, would strongly 
support the fossil disk hypothesis. Until now, the poor spatial resolution of X-ray telescopes have 
made precise positioning of AXPs difficult at best. However, the recently deployed Chandra obser- 
vatory and NewtonXMM mission, with their much improved positioning abilities, should remedy 
this problem and help facilitate optical identifications. In fact, the first possible optical counterpart 
to an AXP has recently been identified (Hulleman, Van Kerkwijk, & Kulkarni 2000). The proposed 
counterpart to 4U 0142 + 61 has a measured flux ratio, Lopt/L^ '-^ 10^^, an order of magnitude 
smaller than what has been predicted for fossil disk models (Perna, Hernquist, &: Narayan 2000). 
Although this measurement seems to, initially, lend more credence to the magnetar hypothesis, it 
should be noted that disk models can not be definitively ruled out either as the emission model of 
Perna, Hernquist, & Narayan (2000) sensitively depends on some poorly understood parameters, 
such as disk size. 



REFERENCES 

Alpar, M.A. 2001, ApJ, 554, 1245 

Angelini, L., Stella, L., & White, N.E. 1991, ApJ, 371, 332 
Baring, M.G. & Harding, A.K. 1998, ApJ, 507, L55 
Bath, G.T. k Pringle, J.E. 1981, MNRAS, 194, 967 
Bondi, H. 1952, MNRAS, 112, 195 
Brown, G.E. 1995, ApJ, 440, 270 

Camilo, F., Kaspi, V.M., Lyne, A.G., Manchester, R.N., Bell, J.F., D'Amico, N.D., McKay, N.P.F., 
& Crawford, F. 2000, ApJ541, 367 

Cannizzo, J.K., Lee, H.M., k Goodman, J. 1990, ApJ, 351, 38 

Chakrabarty, D. 1998, ApJ, 492, 342 

Chatterjee, P., Hernquist, L., k Narayan, R. 2000, (CHN) ApJ, 534, 373 
Chevalier, R.A. 1989, ApJ, 346, 847 
Chevalier, R. 1993, ApJ, 411, L33 

Colpi, M., Geppert, U., k Page, D. 2000, ApJ, 529, L29 



-18- 



Duncan, R.C. 2002 in Soft Gamma Repeaters and Anomalous X-ray Pulsars: The Rome 2000 
Mini- workshop, ed. M. Feroci, & S. Mereghetti, in press (astro-ph/0106041) 

Duncan, R.C, & Thompson, C. 1992, ApJ, 392, L9 

Fabian, A.C. 1975, MNRAS, 173, 161 

FrancischelH, G.J., Wijers, R.A.M.J., & Brown G.E. 2002, ApJ, 565, 471 

Frank, J., King, A. & Raine, D. 1992, Accretion Power in Astrophysics (2d ed.;Cambridge: Cam- 
bridge University Press) 

Ghosh, P., Angelini, L., & White, N.E. 1997, ApJ, 478, 713 

Ghosh, P. k Lamb, F.K. 1979, ApJ, 234, 296 

Hulleman, F., Van Kerkwijk, M.H., & Kulkarni, S.R. 2000, Nature, 408, 689 
Illarionov, A.F. & Sunyaev, R.A. 1975, A&A, 39, 185 
Israel, G.L., Mereghetti, S., & Stella, L. 1994, ApJ, 433, L25 
Kahabka, P. & Li, X.-D. 1999, A&A, 345, 117 

Kaspi, V.M. 2000 in LA.U. Colloq. 177, Pulsar Astronomy - 2000 and Beyond, ed. M. Kramer, N. 
Wex, &, R. Wiebelinski (San Francisco: ASP), 485 

Kudritzki, R.-P. & Puis, J. 2000, ARA&A, 38, 613 

Kunz et al. 1993, A&A, 268, 116 

Leong, C., Kellogg, E., Gursky, H., Tannenbaum, H., & Giacconi, R. 1971, ApJ, 170, L67 
Levine, A., Rappaport, S., Deeter, J.E., Boynton, P.E., & Nagase, F. 1993, ApJ, 410, 328 
Li, X.-D. & Van den Heuvel, E.P.J. 1997, A&A, 321, L25 
Liller, W. 1973, ApJ, 184, L37 

Lin, D.N.C., Woosley, S.E., & Boddenheimer, P.H. 1991, Nature, 353, 827 

Lucke, R., Yentis, D., Friedman, H., Fritz, G., & Shulman, S. 1976, ApJ, 206, L25 

Marsden, D., Lingenfelter R.E., Rothschild, R.E., & Higdon, J.C. 2001, ApJ, 550, 397 

Menou, K., Esin, A.A., Narayan, R., Garcia, M.R., Lasota, J.-P., & McClintock, J.E. 1999, ApJ, 
520, 276 

Mereghetti, S., Chiarlone, L., Israel, G.L., & Stella, L. 2002, in Proc. of the 270 W.E. Heraeus 
Seminar, ed. W. Becker, H. Lesch, &; J. Triimper, in press (astro-ph/0205122) 



-19- 



Mereghetti, S., Israel, G.L., & Stella, L. 1998, MNRAS, 296, 689 

Mereghetti, S., k Stella, L. 1995, ApJ, 442, L17 

Narayan, R. & Yi, I. 1995, ApJ, 452, 710 

Perna, R., Hernquist, L., & Narayan, R. 2000, ApJ, 541, 344 

Price, R.E., Groves, D.J., Rodrigues, R.M., Seward, F.D., Swift, CD., & Toor, A. 1971, ApJ, 168, 
L7 

Primini, F., Rappaport, S., & Joss, P.C. 1977, ApJ, 217, 543 

Pringle, J.E., & Rees, M.J. 1972, A&A, 21,1 

Savonije, G.J. 1979, A&A, 79, 352 

Shakura, N.I. & Sunyacv, R.A. 1973, A&A, 24, 337 

Schreier, E., Giacconi, R., Gursky, H., Kellogg, E., & Tannenbaum, H. 1972, ApJ, 178, L71 

Schwarzchild, M. 1958, Structure and Evolution of the Stars (Princeton:Princeton University Press) 

Thompson, C. 2002 in Soft Gamma Repeaters and Anomalous X-ray Pulsars: The Rome 2000 
Mini- workshop, ed. M. Feroci, Sz S. Mereghetti, in press (astro-ph/0201068) 

Thompson, C., k Duncan, R.C. 1995, MNRAS, 275, 255 

Thompson, C., & Duncan, R.C. 1996, ApJ, 473, 322 

Tjemkes, S.A., Zuiderwijk, E., & Van Paradijs 1986, A&A, 154, 77 

Van Paradijs, J., Taam, R.E., k Van den Heuvel, E.P.J. 1995, AkA, 299, L41 

Van Paradijs, J. k Zuiderwijk, E. 1977, A&A, 61, L19 

Verbunt, F., Wijers, R.A.M.J., k Burm, H.M.G. 1990, A&A, 234, 195 

Wang, Y.-M. 1987, A&A, 183, 257 

Webster, B.L., Martin, W.L., Feast, M.W., k Andrews, P.J. 1972, Nature, 240, 183 



This preprint was prepared with the AAS IM^^X macros v5.0. 



-20- 



TTTI 1 1 — I I I I I l| 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 




I I 1 1 I I I I I I I I I I I I I I I I I I I I I I 

1 10 100 1000 10* 

TIME [Years] 

Fig. 1. — Time evolution of the pulsar spin period as a function of magnetic field strength for a neutron star 
accreting from a fossil disk. Assumed initial conditions are Po = 15 ms, Mo = O.OO6M0, and T = 7/6 (see text). The 
CHN "efficient" propeller torque, given by equation (6), is assumed. We find that a magnetic field strength in the 
range 4 — 10 x 10^^ G is needed for the neutron star to be observed as an AXP. For weaker field strengths, the 
neutron star can not spin down efficiently and is observed as a radio pulsar. When B12 > 3.9, a brief emitter phase is 
followed by an extensive propeller period. The neutron star is observed as an AXP during the "tracking" phase of its 
evolution, i.e. when ~ 1. In our model (as in CHN) the neutron star never enters an extended accretion period. 
The AXP phase is assumed to exist, for the above initial conditions from the transtion time between propeller and 
tracking phases, itrack ~ 10* yr, to the onset of ADAF, 2tADAF ~ 38, 000 yr. 
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Fig. 2. — Final spin period of a neutron star embedded in a fossil disk as a function of opacity parameter, r(/v) for 
various magnetic field strengths. An efficient propeller torque is given by equation (6) and the ADAF transition time, 
tADAP ~ tADAF(r) IS given by equation (5). Five period-opacity curves are shown for neutron star fields ranging 
from 2.5 x 10^^ G to 4.0 x 10^^ G. All curves show the same general trend and the existence of a cutoff opacity 
parameter, Fc. For a given field strength, if F > Fc, the neutron star can not sustain an extended propeller cycle and 
will, consequently, end its life as a radio pulsar. If electron scattering (Fos ~ 1.1875) dominates in the disk, then the 
neutron star can be observed as an AXP provided -B12 > 6.6. For a Kramers disk opacity Fxr = 1.25 S> Fc for all 
B12 < 30. To observe an AXP in a Kramers disk, the minimum field strength necessary is ~ 3.7 X 10^^ Gauss. In 
their analytic model, Chatterjee, Hernquist, & Narayan (2000) assumed F = 7/6 = 1.167. 




Fig. 3. — Spin evolution for a neutron star embedded in a fossil disk for three different propeller torques. Initial 
conditions for the various torques are otherwise fixed:Mo = O.OO6M0, Pq = 15 ms, Bo = 7.5 x 10^^ Gauss, and 
T(k) = 7/6. Employing a CHN propeller torque yields, after a sharp propeller phase, an AXP in the tracking phase 
near t ~ 1.5 x 10* yr. The neutron star remains at quasi-cquilibrium until 2tADAF ~ 3.8 x 10'* yr with a final spin 
period, 13.5 s. For both the energy (E) and angular momentum (J) propellers, the neutron star never leaves the 
propeller phase of its evolution and the final spin period is ~ 0.2 s. In both cases, the neutron star is observable as a 
radio pulsar. The early evolutionary deviations for the energy and angular momentum propellers are due to the lack 
of an emitter phase for the former. At later times, their spin evolution is consistent. 
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Fig. 4. — Spin evolution models for SMC X-1 assuming various propeller torque models. The progenitor system, 
assumed to be an 09 V star and a newly born neutron star, have the following initial properties: Mi = 20Mq, 
Ri = WRq, B„ = 10^'^ G, Po = 15 ms. Mass is transfered by a radiation-driven wind during the (non-conservative) 
main sequence phase of the system's evolution. At t = 5.6 Myr, Hydrogen shell burning initiates in the primary, 
instigating atmospheric Roche-lobe overflow. If an energy propeller is assumed, the calculated spin period during 
the suporgiant phase ranges from 5.3 s, at the onset of the heavy accretion cycle, to 0.99 s at t = 6.2 Myr. For the 
angular momentum propeller, P{ ranges from 2.15 s to 0.52 s. Assuming the efficient propeller torque of Chatterjee, 
Hernquist, & Narayan (2000) yields a neutron star with a spin period ranging from 15.9 s to ~ 100 s. The observed 
spin period of SMC X-1 = 0.71 s is illustrated for reference. 



